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Abstract 

As  part  of  the  Galileo  project  of  Early  Trials  on  Time  Synchronization  Techniques  and 
Calibration  Issues,  the  (near)  zero-baseline  GPS  Common-View  Time  Transfer  and  Direct 
Measurement  (using  a  Time  Interval  Counter)  results  between  the  time  scales  of  two  co-located 
clocks  during  2001  were  analyzed.  This  was  done  to  identify  the  obtainable  uncertainty  levels 
with  these  techniques.  The  result  of  the  analysis  is  presented.  The  baseline  of  the  antennas  of 
the  two  single-channel  NBS-type  GPS  C-V  receivers  was  2  m,  so  multi-path  influence  is  not 
excluded.  The  uncertainties  from  statistical  origin  (Type  A  evaluation)  were  determined  by 
means  of  MDEV  and  TDEV.  The  uncertainties  from  other  sources  (Type  B  evaluation)  were 
determined  from  calibrations.  The  TDEV  of  the  GPS  C-V  result  was  less  than  0.4  ns  at 
averaging  times  >  1  day,  and  at  shorter  averaging  times  TDEV  was  about  1  ns. 


1  INTRODUCTION 

The  work  described  in  this  paper  was  done  in  the  framework  of  the  Galileo  project  as  part  of  an  ESA 
Early  Trials  study  on  Time  Synchronization  Techniques  and  Calibration  Issues.  The  goal  was  to  verify 
the  1  ns  GPS  Common-View  (C-V)  clock  synchronization  accuracy  assumption  in  a  laboratory  (“best 
case”)  co-located  clock  comparison  experiment. 

For  this  purpose  during  2001  a  laboratory  test  has  been  performed  at  VSE  to  compare  the  time  scales 
from  two  co-located  cesium  clocks  by  direct  measurement  using  a  time  interval  counter  (TIC)  and  by 
means  of  two  co-located  (near  zero  baseline)  GPS  C-V  receivers.  The  latter  is  the  “best  case”  situation  of 
GPS  C-V,  because  the  high  correlation,  and  hence  first-order  cancellation,  for  environmental  effects  such 
as  temperature,  humidity,  troposphere,  ionosphere,  rainfall,  etc.  The  difference  of  the  results  of  the  two 
methods  is  given  and  discussed. 


2  GENERAL 

2,1  Clocks  and  Time  Scales  Used 

NMl-VSE  uses  a  commercial  high-performance  cesium  clock  (HP5171A)  followed  by  a  phase  micro¬ 
stepper  for  UTC  (VSE),  which  serves  as  the  direct  reference  for  the  first  GPS  receiver  (VSEOl).  A 
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second  high-performance  cesium  clock  (HP5171A)  followed  by  a  microstepper  generates  a  backup  time 
scale  called  UTC  (VSB)  and  serves  as  reference  for  the  second  GPS  receiver  (VSL18). 

2,2  Accuracy  and  Uncertainty  Expression 

The  ISO  Guide  to  the  Expression  of  Uncertainty  in  Measurement  (GUM)  [4]  and  the  EAL  document 
EAE-R2  [5]  discern  two  types  of  uncertainties  related  to  measurements:  the  uncertainty  evaluated  from 
the  statistical  analysis  of  the  data,  denoted  Type  A  evaluation,  and  indicated  by  Ua;  and  the  uncertainty 
evaluated  from  other  causes,  also  denoted  Type  B  evaluation,  indicated  by  Ub.  Type  A  evaluation 
includes  contributions  from  random  nature.  Type  A  in  the  field  of  Time  and  Frequency  is  mostly 
expressed  as  Allan  deviation  (ADEV),  Modified  Allan  deviation  (MDEV),  and  Time  deviation  (TDEV) 
[1,  2].  Type  B  evaluation  includes  estimates  of  contributions  from  systematic  nature,  derived  from 
calibrations,  closing  errors,  etc. 

The  total  combined  uncertainty  indicated  by  Uc  is  calculated  as  the  root  of  the  sum  of  the  quadratics  of  Ua 
and  Ub.  To  obtain  a  confidence  level  of  95%,  a  coverage  factor  k  =  2  is  generally  used.  In  the  form  of  a 
formula:  Uc  =  2  x  (V(V(  Ua)+(V  (  Ub))).  The  term  “accuracy”  should  be  avoided  for  measurement  results. 


3  DIRECT  CLOCK  MEASUREMENT  (SHORT  CABLE  TIC) 

3.1  Description  of  the  Direct  Clock  Measurement  Setup 

The  differences  between  the  two  time  scales  VSE  and  VSB  are  measured  once  each  hour  directly  with  a 
time  interval  counter  (TIC).  A  TIC  is  a  device  which  measures  the  time  interval  between  the  signal 
applied  to  its  start  input  compared  to  the  signal  applied  to  its  stop  input  (see  Figure  1).  The  resolution  of 
the  measured  and  displayed  time  interval  in  the  best  TIC  can  be  10  picoseconds  (see  Figure  2),  while  the 
Type  B  uncertainty  (accuracy)  can  be  50  ps  to  100  ps.  The  time  scale  signals  and  the  time  outputs  from 
the  used  GPS  receivers  are  steep  pulses  occurring  once  per  second  (1  pulse  per  second,  1  PPS).  The  TIC 
and  the  associated  coaxial  switching  system  with  cables  were  calibrated.  The  TIC  was  calibrated  for 
internal  start-stop  delay  by  connecting  the  start  signal  to  both  inputs.  The  delay  in  the  remaining  system 
was  determined  by  connecting  the  start  signal  to  the  input  of  the  system,  where  the  UTC  (VSB)  signal 
was  connected.  The  total  correction  on  the  direct  measurement  data  (for  the  VSB  signal)  to  be  applied  is 
-34.6  ns.  Outliers  in  the  direct  measurement  >  200  ns  have  been  removed,  except  for  known  steps  of  a 
multiple  of  100  ns.  These  steps  have  been  corrected. 

3.2  Corrected  Direct  Data 

The  direct  VSE-VSB  data  are  corrected  for  known  time  steps,  and  the  result  is  shown  in  Figure  3.  The 
last  part  of  the  Figure  3  shows  a  straight  line.  This  is  not  an  error,  but  during  this  period  the  second  clock 
B  was  been  replaced  by  Clock  A,  so  during  this  period  a  common  clock  was  used  for  the  direct  clock 
measurements. 

3.3  Characterization  and  Analysis 

The  MDEV  and  TDEV  of  the  direct  clock  difference  results  are  shown  in  Figure  4.  The  MDEV  shows  a 
straight  line  with  a  slope  of -A  dependence  on  x  from  3600  s  (1  h)  to  about  4-10^^  s  (4  days).  Then  some 
flattening  occurs. 
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For  TDEV  it  is  also  a  straight  line,  but  with  a  slope  of  +V2  dependency  on  x  from  200  ps  at  1  h  to  2.3  ns  at 
4  days.  Please  realize  here  that  these  results  are  not  giving  only  the  time  transfer  properties;  they  include 
also  the  properties  of  the  two  clocks. 


4  THE  CO-LOCATED  GPS  C-V  METHOD 

4.1  Indirect  Measurement  Setup  Using  GPS  Receivers  in  Common-View 

Using  the  GPS  Common-View  technique,  two  UTC  time  scales  have  been  compared  using  two  GPS  time 
transfer  receivers.  This  co-located  GPS  (Galileo)  Common-View  time  transfer  method  (Figure  1)  is 
identical  to  the  general  GPS  C-V  method.  The  difference  is  that  the  receivers  and  their  antennas  are  now 
co-located  within  2  m,  thus  near  zero-baseline.  In  such  case  the  ionosphere  and  troposphere  delays  are 
identical  and  cancel  completely  in  the  GPS  C-V  result.  The  only  differences  are  related  to  the  multi-path, 
antenna,  cables,  and  receiver  delay  differences.  So  this  method  can  be  used  to  measure  the  optimum,  but 
realistic  time  comparison  properties  (realistic  upper  performance  limit)  for  GPS  C-V  using  these  two  co¬ 
located,  mutually  independent  receivers. 

The  two  GPS  receivers  VSLOl  and  VSL18  used  for  this  experiment  at  VSL  are  Time  Transfer  receivers 
of  the  NBS-type,  single-channel  C/A  code  tracking. 

4.2  Co-LOCATED  GPS  C-V  Data  Results 

The  1-year  raw  data  from  VSLOl  and  from  VSL18  have  been  corrected  for  known  jumps.  Each 
individual  data  point  of  the  13-minute  sessions  according  to  the  BIPM  schedule  for  both  receivers  is 
shown  in  Figure  5. 

The  clock  difference  data  of  the  VSL-GPS  and  VSB-GPS  data,  named  (VSL-VSB)gps  data,  are  also  shown 
in  Figure  5  (green,  upper  trace).  The  last  part  of  this  trace  shows  a  straight  line.  This  is  not  an  error,  but 
during  this  period  after  MJD  52225  the  second  clock  B  generating  time  scale  VSB  was  been  replaced  by 
the  Clock  A  that  generated  also  the  first  time  scale  UTC  (VSL),  so  a  common  clock  was  used  for  both 
GPS  receivers  during  this  period.  The  cesium-beam  tube  of  Clock  B  was  at  its  end-of-life. 

4.3  Characterization  and  Analysis 

The  MDEV  and  TDEV  of  the  clock  difference  data  are  shown  in  Figure  6.  The  TDEV  of  this  near  zero- 
baseline  GPS  C-V  time  transfer  result  slowly  increases  with  averaging  time  x  from  1  ns  at  1  h  to  2  ns  at  2 
days  and  reaches  1 0  ns  at  x  of  40  d. 

The  frequency  transfer  result  between  the  clocks  shows  a  modified  Allan  deviation  of  2T0'’"'  at  1  day  and 
has  a  floor  of  5-10'’^  at  >  20  d.  Please  realize  that  these  results  are  not  giving  only  the  time  transfer 
properties;  they  include  also  the  properties  of  the  two  clocks. 
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5  DIFFERENCE  GPS  C-V  AND  DIRECT 

5.1  Interpolated  and  Corrected  Differential  Data 

The  difference  between  the  direct  method  and  GPS  C-V  is  shown  in  Figure  7.  This  difference  is  ((VSL- 
VSB)  direct  -  (VSL-VSB)  GPS.  After  MJD  52225,  the  VSL  clock  was  replaced  the  VSB  clock,  so  from 
that  time  on  it  is  common  clock  as  can  seen  from  the  straight  noiseless  (VSL-VSB)  direct  difference.  The 
yellow,  lower  trace  is  the  “double”  difference.  There  is  no  difference  visible  during  the  period  of 
common  clock,  so  both  clocks  have  similar  stabilities. 

5.2  Characterization  and  Analyses 

The  result  for  the  difference  of  the  direct  minus  the  GPS  C-V  results  are  shown  in  Figure  8.  It  can  be 
seen  that  the  TDEV  is  at  1  ns  level  for  x  <  40  000  (0.5  d),  then  becomes  for  x  >  1  day  rather  constant  at 
the  400  ps  level  (ua)  for  x  up  to  20  d.  The  latter  reflects  mainly  the  differential  stability,  the  sum  of  all 
other  long-time  uncorrelated  factors,  of  the  two  GPS  receiver  results,  while  the  first  shows  multi-path, 
diurnal  differential  temperature  dependency,  and  other  short-time  uncorrelated  sources  such  as  basic 
receiver  phase  noise. 


6  CONCLUSIONS 

6.1  Type  A  Uncertainty  Evaluation 

For  near  zero-baseline  GPS  C-V,  the  uncertainty  from  statistical  origin  Type  A  (ua)  is  for  k  =  1  equal  to 
the  TDEV,  and  so  is  1  ns  for  x  <  40  000  s  (=  0.5  d),  and  400  ps  for  x  >  1  to  20  days. 

6.2  Type  B  Uncertainty  Evaluation 

The  average  difference  between  the  direct  method  and  the  GPS  C-V  method  for  time  transfer  without  any 
correction  was  +  46  ns.  Application  of  the  necessary  correction  of  -34.6  ns  for  the  VSB  direct  data  brings 
the  difference  (direct  -  GPS  C-V)  to  the  value  46  -  34.6  =  +1 1.4  ns. 

Application  of  a  correction  of +1 1.7  ns  on  the  BIPM  calibration  of  VSEOl  in  1998  [3]  to  the  VSEOl  data 
and  a  correction  of +1.8  ns  to  the  VSL18  data  give  a  correction  for  the  GPS  C-V  results  of +1 1.7  -1.8  = 
+9.9  ns.  This  reduces  the  difference  (direct  -  GPS  C-V),  the  closing  error,  to  +1 1.4  -  9.9  =  +1.5  ns  (=  Ub). 
This  is  an  excellent  result  for  Type  B  evaluation  Ub. 

6.3  Combined  Uncertainty  Evaluation 

Taking  the  direct  measurement  as  reference,  then  the  total  combined  uncertainty  Uc  for  the  co-located, 
near-zero-baseline  GPS  C-V  can  be  estimated  to  be  for  k  =  2,  at  a  95%  confidence  interval  (see  Section 
2.2),  up  to  0.5  d  averaging  time:  3.6  ns,  and  for  1  to  20  day  averaging  time:  3.1  ns. 
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6,4  Remarks 

Improvements  in  the  GPS  C-V  results  could  result  from  application  of  postprocessed  data:  satellite  orbits, 
measured  ionosphere  data,  measured  troposphere  data,  and  also  from  improvement  of  temperature 
sensitivity  and  multi-path  in  receiver  equipment  and  cables. 
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Co-located  Clocks  through  1:GPS  Common-View 
and  2:  Direct  Time  Interval  Counter  measurement 


Time  difference  1  (A-B)  =  (A-C)  -  (B-C)  +  delay  diff. 
Time  difference  2  (A-B)  =  TIC  reading 


J 


Figure  1.  Co-located  measurement  setup. 


Stability  of  Time  Interval  Counter  delay  measurement 


tau  (s) 


Figure  2.  Time  interval  counter  MDEV  and  TDEV. 
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Clock  Difference  (VSL-VSB)  by  Direct  Measurement 


I  ■  (VSL-VSB)direct  | 


Clock  Difference  (VSL-VSB)  by  Direct  Measurement 


I - MDEV(direct) - TDEV(direct)  | 


Figure  4.  Characterization  of  VSL  -  VSB  clocks  by  direct  measurement. 
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Clock  Difference  (VSL-VSB)  by  co-located  GPS  C-V 


I  »  VSL-GPS  ■  VSB-GPS  ■  (VSL-VSB)gps  | 
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Figure  5. 


GPS  C-V  data  and  the  clock  difference  (VSL  -  VSB)  by  GPS  C-V. 


Clock  difference  (VSL-VSB)  by  co-located  GPS  C-V 


tau  (s) 

Figure  6.  Characterization  of  clock  difference  (VSL  -  VSB)  using  co-located  GPS  C-V. 
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Clock  Difference  (VSL-VSB)  by  GPS  C-V  and  Direct  Measurement 

I  •  (VSL-VSB)g^s  ^  (VSL-VSB)direct  (VSL-VSB)direct-(VSL-VSB)gps  I 


Difference  ((VSL-VSB)direct  -  (VSL-VSB)co-located  GPS  C-V) 


1^— MDEV  ((VSL-VSB)direct-(VSL-VSB)gps)  TDEV  ((VSL-VSB)direct-(VSL-VSB)gp5)  | 


tau  (s) 

Figure  8.  Characterization  of  the  difference  between  co-located  GPS  C-V  and  direct  measurements. 
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Questions  and  Answers 

CHRISTINE  HACKMAN  (University  of  Colorado):  My  question  is  on  your  time  deviation  plot,  after 
about  a  day,  it  looks  like  you  have  flicker  PM  at  about  400  ps.  I  was  wondering  if  you  know  what  the 
source  of  that  is. 

GERRIT  de  JONG:  I  think  these  are  the  remaining  differences  for  the  long  term.  So  that  could  be  the 
temperature  coefficients  of  both  systems.  It  could  be  in  first  order  the  same,  but  not  exactly  the  same.  So 
then  you  could  have  that  remaining  fluctuation,  depending  on  temperature  or  any  other  thing,  which  does 
not  compensate  fully  in  the  two  receivers.  They  are  of  the  same  kind,  the  same  design,  but,  of  course, 
there  are  individual  differences  in  the  antennas,  mainly,  because  there  you  have  the  biggest  environmental 
influence  that  could  still  be  present.  One  antenna  is  directly  coming  from  NIST  at  that  time.  The  other  is 
an  Allen  Osborne  antenna.  So  they  are  not  from  the  same  unit,  and  there  could  be  some  remaining 
difference. 
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